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We propose a few selected experimental approaches to show that new-generation instruments can 
give a direct access to significant observables on the density dependence of the symmetry energy in 
the nuclear equation of state. The form of such dependence is investigated within the Stochastic 
Mean Field model, coupled to a secondary-decay treatment. 



I. INTRODUCTION 

The exotic-beam facilities which already exist, or 
which are under construction, invite to focus on a long- 
standing experimental challenge: heavy-ion collisions 
induced by exotic nuclei at low and intermediate en- 
ergies will be a probe for the properties of isospin- 
asymmetric nuclear matter In particular, numer- 
ous next-generation experiments will be dedicated to the 
study of the density behaviour of the symmetry energy. 

Such quantity defines the isospin-asymmetric part of 
the equation of state (EOS), which can be deduced from 
different phenomenological parametrisations [2- 4]. The 
corresponding predictions can be grouped schematically 
into two main forms: either a 'stiff' or a 'soft' dependence 
of the symmetry energy as a function of the nuclear den- 
sity d, [f|. At present, considerable effort is dedicated to 
pin down the density dependence in the regime of rela- 
tivistic heavy- ion reactions [X, H, @ , where supersatura- 
tion densities are accessible. Such conditions determine a 
significant difference between the two forms for the den- 
sity dependence but the uncertainties in the hadron ef- 
fective mass splittings and momentum dependence in the 
isovector channel introduce new degrees of freedom, en- 
hancing the complexity of the problem Q. 

These complications do not arise in the Fermi-energy 
regime, where the functional form of the symmetry en- 
ergy at subsaturation densities should in principle be ac- 
cessible. Furthermore, these conditions give access to 
the interesting issue of the influence of cluster correla- 
tions [&, 9f. At Fermi energies, the difficulty is that the 
difference between the theoretical predictions for the den- 
sity dependence is less pronounced than at the large su- 
persaturation densities which can be probed at rclativis- 
tic energies: in particular, it is not always clear whether 
an effective discrimination between different degrees of 
asy-stiffness is really possible. On the one hand, this 
is due to the dependence of the predictions on the ex- 
isting transport models 0, 0, [H|. On the other hand, 
besides the detection limitations, also the secondary de- 
cay contributes in deforming EOS-scnsitivc isotopic ob- 



servables [Ill4l3j : we will discuss how this can affect the 
comparison protocol between model and experiment. 

In this respect, we need to identify isotopic observ- 
ables which, first of all, manifest a significant sensitiv- 
ity to the change between the different asy-EOS forms 
and, in addition, constitute robust experimental observ- 
ables against the effect of secondary decay and detection 
limitations [HI EMI, is necessary. We also indicate 
that, recently, such study has profited from the results 
of isospin-transport experiments, which already imposed 
theoretical constraints [fl [l7l - l2~l| . This work is aimed 
to give a schematic guideline for measuring a selection of 
isospin observables with new-generation isospin-sensitive 
instruments, in comparison with former isospin-blind de- 
vices. In particular, within our simulation protocol, we 
focus on a forthcoming Air detector, FAZIA (22j, explic- 
itly planned for measuring isospin observables and its 
isospin-blind ancestor, INDRA [23|,[24|- We only focus on 
these two detectors because they are well suited for rep- 
resenting two opposite situations. However, it is evident 
that they can not resume the several experimental strate- 
gies which exist or are planned and which are based on 
different innovative tools for accessing the isospin observ- 
ables. In particular, the An detector CHIMERA (2f| is a 
new-generation detector which is rapidly improving and 
evolving between these two extremes and its experimen- 
tal results already mark the way towards future isospin- 
sensitive devices. 



II. DESCRIPTION OF THE MODEL 

In this paper, a collision system is described within 
the Stochastic Mean Field (SMF) model. This is a 
time-dependent semi-classical mean-field model where 
nucleon-nucleon collisions as well as fluctuations are 
taken into account in the Boltzmann-Uehling-Uhlenbeck 
formalism (see [26| and references therein); the refer- 
ences [1, [27} give theoretical and numerical details. Such 
model is reliable in the Fermi-energy regime, for an in- 
cident energy ranging between about 10 and 200 MeV 
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FIG. 1: (Color online) Left. Asy-soft and asy-stiff forms for 
the symmetry energy. Right. Corresponding derivative with 
respect to the density in the low-density region 



per nucleon. In particular, due to the realistic isospin- 
dependent mean-field and the introduction of density 
fluctuations, the model gives a reliable description of 
heavy and intermediate-mass fragment properties; since 
elastic collision processes are accounted for, the model 
also describes the production of prccquilibrium neutrons 
and protons reliably; nevertheless, the description of light 
particles demands correlations which are out of the scope 
of the model. 

In the calculations, we adopt a parametrisation which 
gives the same properties as the SkM*[28j force for 
nuclear matter: these are a saturation density po = 
0.16 fm -3 and an incompressibility modulus = 
200 MeV which corresponds to a soft equation of state. 
We adopt two different prescriptions for the behaviour of 
the symmetry energy S, which arc respectively: 
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with the parameters a = 12.7 MeV, b = 36McV, c = 
481.7 MeV fm 3 , and d = -1638.2 MeV fm 6 . As illus- 
trated in fig. [TJ the first is an 'asy-stiff' form, for which 
the potential symmetry term linearly increases with nu- 
clear density, while the second is an 'asy-soft' form, corre- 
sponding to a flatter behaviour of the potential symmetry 
energy around and below normal density. More details 
are provided in Ref.Q. 

Isospin effects originate from the fact that neutrons 
and protons experience different forces. In particular, 
the difference between the neutron and proton currents 
in — jp, that develop in presence of asymmetry (VI) 
and/or density (V p) gradients, is stricly connected to the 
strength of the symmetry energy (and of its derivative). 
In fact, within a simple hydrodynamical picture, one can 
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Hence, in presence of asymmetry gradients (diffusive pro- 
cesses) we test essentially the strength of the symme- 
try energy while, when density gradients are encoun- 
tered along the dynamical path, we observe "isospin mi- 
gration" towards the low density regions, ruled by the 
derivative of the symmetry energy. 

Within the SMF model, the systems 68 Ni+ 68 Ni, 
58 Ni+ 68 Ni and 58 Ni+ 58 M at 15 and 4Q MeV per nucleon 

were simulated in the present work for a continuous dis- 
tribution of impact parameters b, which evolves as bdb. 

Fig. [2] shows the reaction mechanisms observed as 
a function of the reduced impact parameter 6 re( j (nor- 
malised to the sum of the nuclear radii of the target 
and projectile nuclei). Fusion reactions (F) were iden- 
tified from the presence at the asymptotic time of a 
unique fragment having bigger mass than the projectile 
(or the target, equivalently) ; the quasi-elastic channel is 
defined as non-fusion events where the sum of the charges 
of the two largest fragments exceeds a large fraction of 
the charge of the projectile (or the target, equivalently; 
this fraction is adjusted to 90% and 85%, for the inci- 
dent energies of 15 and 40 AMeV, respectively); the neck 
contribution (neck) is recognised from events which do 
not belong to the previous categories and where at least 
three fragments have larger charge than helium; all the 
other events are interpreted as dissipativc binary reac- 
tions (DI). This figure shows that a centrality selection 
corresponds only approximately to a selection of the re- 
action mechanism. This is especially true for the neck 
events, which never dominate the total cross section. 
Since the experimental apparatus does not deliver the 
same response to the different reaction mechanisms, this 
may result in an uncontrolled bias if unfiltcred simula- 
tions are confronted to the experimental data. 
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FIG. 2: (Color online) Probability for different reaction 
mechanisms as a function of the reduced impact parameter for 
the systems 68 Ni+ 68 Ni at 15 (left panel) and 40 (right panel) 
AMeV. The asy-soft (here used) and asy-stiff (not shown) 
forms give very close results. 
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A. Evolution of the isospin content and effect of 
the secondary decay 

In the calculations, the dynamical stage is followed 
till the time t = 260 fm/c. In order to deal with the 
full range of impact parameters which are simulated, the 
fragment properties are evaluated as soon as the system 
breaks up into pieces within this interval of time, in each 
event. This procedure implies that events with only one 
fragment (incomplete fusion) are followed until the final 
(longer) time, leading to a reduction of the calculated 
size and excitation energy of the compound system, due 
to nucleon emission. The output of the SMF calculation 
is then coupled to the evaporation model GEMINI [29[ 
(the 2003 release is used). This coupling is necessary be- 
cause of two reasons. First, secondary decay extends over 
a much longer time than the numerical stable interval 
of time of any transport code; Second, even molecular- 
dynamics models, which in principle can describe the pro- 
duction of complex light particles, can not reproduce the 
statistical branching ratios of compound-nucleus emis- 
sion m. 

We shall emphasise that in the absence of secondary 
decay, a basic quantity like the isotopic composition 
of intermediate-mass fragments (IMF) would be suffi- 
cient to distinguish between the asy-stiff and asy-soft be- 
haviour of the symmetry energy, even in stable systems. 
This is shown in fig. [3] (upper spectra labelled 'SMF'), 
where the average neutron number to atomic number 
ratio of the biggest fragment in all events not issued 
of fusion is plotted as a function of the corresponding 
longitudinal component of the rapidity in the centre of 
mass, relative to the projectile (target) rapidity for posi- 
tive (negative) rapidities, for the two systems 58 Ni+ 58 Ni, 

58 Ni+ 68 Ni and 68 M+ 68 Ni &t lg and 4 q ^ MeV . In the 

symmetric systems, the fragments are more neutron rich 
in the asy-stiff case and the effect enhances for the more 
neutron-rich system. This scenario can be attributed 
to the lower value of the symmetry energy below normal 
density in the asy-stiff case, with respect to the asy-soft 
form (sec fig. [TJ. In the asymmetric system, in addi- 
tion to nucleon emission, isospin diffusion takes place 
through the low density interface between the two re- 
action partners. As illustrated by eq(31 diffusion is a 
process of isospin equilibration between the two asym- 
metric reaction partners which acts more effectively for 
higher values of the symmetry energy; at low densities, 
this process is therefore more effective for a softer form 
of the EOS, which is in this case given by the asy-soft 
parametrisation (see fig. [1}. As a consequence, the pro- 
jectile and target sides in the diagram of fig. [5] are more 
similar in this case, and the difference between the two 
paramctrisations is particularly relevant at midrapidity, 
for the system at 40 A MeV 

If the reaction ended at this stage in the laboratory, 
there would be no need to introduce complex differential 
observables [13], or to exploit exotic radioactive beams 
for the purpose of this investigation. 



Nevertheless, unfortunately, the secondary decay 
makes the discrimination between the two parametrisa- 
tions difficult: in stable and moderately exotic systems 
the secondary decay can wash out most of the discrim- 
inating signals; this effect is illustrated in fig. [3] (lower 
spectra labelled 'SMF+GEMINF), where the hot system 
described within the SMF model is made decay by the 
GEMINI model. 

We can not exclude that such a prominent effect may 
depend on the model or the system chosen; different con- 
tributions of secondary decay have been reported in the 
literature [TlTflil and no systematic study has been 
performed to our knowledge. However, this prologue 
should drive the attention on the fact that the direct mea- 
surement of observables related to the equation of state is 
certainly challenging. In particular, the evaporation pro- 
cess has the effect of modifying and even smearing out 
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FIG. 3: (Color online) Variation of the isotopic composition 
of the biggest fragment (excluding fusion events) as a function 
of the rapidity in the centre of mass relative to the projectile 
(target) rapidity for positive (negative) rapidities calculated 
with both the asy-stiff and asy-soft parametrisations before 
(SMF, upper spectra) and after (SMF+GEMINI, lower spec- 
tra) secondary decay for the systems 58 Ni+ 58 Ni 58 Ni+ 68 Ni 
and 68 Ni+ 68 Ni at 15 AMeV (upper panels) and 40 AMeV 
(lower panels), for all impact parameters. The dots indicate 
fusion nuclei. The dashed lines indicate the average isotopic 
composition of the colliding system. 
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the observables wc wish to investigate. Beside achieving 
an increasingly more precise understanding of the decay 
sequence (especially in exotic systems), we should also 
define efficient strategies in choosing the nuclear system, 
adopt new observables and refine the criteria of event 
selection to reduce the effect of the secondary decay, so 
as to preserve significant signatures. These requirements 
also impose to develop innovative experimental strate- 
gics. 



B. Detection 

In the present report we suggest some possible solu- 
tions which would require the use of a dedicated detec- 
tor device. Such detector, which is under development 
within the FAZIA project ('Four-7r A and Z Identifica- 
tion Array' [22j ), is a 47r-array of telescopes, each one 
designed for measuring the kinetic energy, the nuclear 
charge and the mass of the intercepted fragments (the 
mass measurement over an extended mass range consti- 
tutes the main innovation for such device). An addi- 
tional calculation filters the results produced by the re- 
action model (SMF+GEMINI) by accurately simulating 
the functioning of the FAZIA detector, including the de- 
tailed geometry and the response of the different detec- 
tion modules [31j. In particular, in the case of FAZIA, the 
telescope is composed by two silicon detectors of 300 /iin 
and 500 /im respectively and a Csl detector; if a particle 
is arrested in the second silicon detector, the kinetic en- 
ergy and the time of flight arc measured. In this case, 
the particle can be identified in nuclear charge and, if 
possible, in mass by the correlation of the energy versus 
the energy loss; in the simulation the measured nuclear 
charge is then set equal to the theoretical value, and the 
mass is set equal to the theoretical value up to phosphor 
and, for higher elements, it is deduced from the kinetic 
energy and the time of flight. 

The functioning of the former-generation detector IN- 
DRA [23|, Hil is also simulated for comparison. FAZIA 
differs from INDRA because the granularity is larger 
(with a substantial gain in angular resolution) and be- 
cause the mass, in addition to the nuclear charge, can be 
measured. When necessary, we will also simulate in the 
calculation the existence of neutron detectors imagined as 
a belt of DEMON [32[ telescopes covering all azimuthal 
angles and disposed on the longitudinal plane (containing 
the beam axis) . For such a configuration, we do not take 
into account the effect of the charged-partiele detector on 
the neutron detection in the DEMON modules. 

In this report, we compare two exclusive experimen- 
tal approaches for reaction experiments with exotic nu- 
clei. The first employs INDRA (with, if needed , the 
addition of DEMON modules), where the mass of the 
fragments (and any isotopic observable) is not measured, 
but deduced from the nuclear charge through the use of 
a paramctrization (EPAX, which supposes that the cold 
nuclei are attracted towards the residue corridor) , which 



is closely compatible with the decay model GEMINI and 
which was adjusted to mostly stable (or moderately ex- 
otic) nuclei. The second approach is FAZIA (with, even- 
tually, the addition of DEMON modules), which delivers 
a direct measurement of the masses, independently of the 
exoticity of the system. Hence the comparison with the 
response of INDRA is useful to appreciate up to which 
extent, for the experiment and the observables under ex- 
amination, fragment masses may deviate from the residue 
corridor and, once measured by FAZIA, bring informa- 
tion on the isovector term of the nuclear interaction. 



C. Statistics 

The number of events considered in the transport cal- 
culations are 1000. Each one of these events is then used 
to produce ten decay paths. For the purpose of simulat- 
ing the experiment, since the geometrical efficiency in the 
selection of complete events could impose a severe reduc- 
tion of the statistics, the number of events is then largely 
increased by considering rotations around the beam axis. 
At the end, the simulated values, within the uncertain- 
ties, correspond to one standard deviation of the mean, 
calculated for a statistics of 10000 events (1000 transport 
events multiplied by 10 decay paths) and considering the 
full error propagation. We would like to stress that the 
simulated results we give should be considered as a lower- 
limit expectation with respect to a real experiment. The 
reason is that in a real experiment we can profit from 
advanced tools of event selection and develop the data 
analysis down to fine details, while in the simulation the 
reduced statistics allows to simulate the data analysis 
only up to a limited degree of accuracy. This limita- 
tion also propagates to the error bars of the simulated 
results. 

For this same reason, we did not exploit the full set 
of simulations. In fact, in some cases (especially for the 
neutron-poor systems at the lower beam energy) the un- 
certainties were too large, as compared to the different 
predictions of the two asy-EOS. We stress again that this 
would not correspond well to the experimental situation, 
which is already successfully progressing on the study of 
isospin effects with stable beams. Moreover, for the pur- 
pose of comparison with experimental data, the accuracy 
of the simulated results can be improved by increasing 
the number of SMF events considered. 



III. SIMULATIONS: FROM THE SELECTION 
OF THE IMPACT PARAMETER TO THE 
ISOTOPIC OBSERVABLES 

Experimentally, the impact parameter can be deduced 
from the transverse energy per nucleon E tT of light par- 
ticles with Z < 3 (it should be precised that in the SMF 
simulations no particles with Z = 2 arc produced; they 
are generated in the decay, by GEMINI). The correlation 
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between the impact parameter and the transverse energy 
per nucleon is illustrated in fig.0]for 58 Ni+ 58 Ni at 15 and 
40 AMcV. At 40 AMeV, the correlation extends over all 
impact parameters except for 6 re d < 0.2, where fusion be- 
comes the dominant mechanism (see fig- [2]) . At 15 AMeV, 
the correlation is still valid for peripheral collisions, but 
the sensitivity is gradually lost when entering the fusion 
regime (& re d < 0.55) and the dissipation becomes total 
for complete fusion (around & r cd < 0.3). In fig. [4] the 
measurement of such correlation is simulated for the two 
detector arrays INDRA and FAZIA: the measured trans- 
verse energy per nucleon is divided by the geometric ef- 
ficiency e gco (INDRA: e goo = 0.88, FAZIA: e goo = 0.78). 
Both INDRA and FAZIA perform equivalently well. We 
shall emphasise that this implies that this same conclu- 
sion extends to any assembly where a part of one array is 
replaced by the corresponding part of the other, without 
reducing the whole angular coverage. 

Profiting of this selection, we focus on the study of 
the isotopic composition of the largest fragments. This 
study, presented in fig. [5j is one example intended to 
answer the question whether the advantage of disposing 
of exotic systems could counterbalance the effect of the 
evaporation process. In fig. [SJ three systems 58 Ni+ 58 Ni, 
58 Ni+ 68 Ni and 68 Ni+ 68 Ni at two incident energies are 
studied; to keep an approximate link with the underly- 
ing reaction mechanism (see fig. [5]), central events are 
defined as 6 rc d < 0.2 (6 ro d < 0.55) for the incident en- 
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FIG. 4: (Color online) Correlation between the impact pa- 
rameter and the transverse energy per nucleon for 58 Ni+ 58 Ni 
at 15 and 40 AMeV, as simulated by SMF+GEMINI (solid 
lines). For comparison, the same correlation is simulated for 
the two experimental approaches INDRA and FAZIA, where 
the transverse energy per nucleon is further divided by the 
geometric efficiency e gco . The interval in impact parameter 
which corresponds to fusion is indicated. The statistical un- 
certainties of the calculation SMF+GEMINI (asy-soft form) 
are indicated by the coloured bands; they are reflected in the 
filtered data with comparable magnitude (Not indicated in 
the figure, for better visibility) 
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FIG. 5: (Color online) Isotopic composition of the largest 
fragments as a function of the isotopic composition of the 
system at 15 (upper panels) and 40 AMeV (lower panels), 
under the constraint of selecting either central (left panels) or 
mid-peripheral events (right panels). This observable is sim- 
ulated for the hot system (SMF - asy soft) and for the cold 
system (SMF+GEMINI). For comparison, the same correla- 
tion is simulated for the two experimental approaches INDRA 
and FAZIA. 



ergy of 40 (15) AMeV; mid-peripheral events are defined 
as 0.25 < 6 ro d < 0.6 (0.6 < 6 rcd < 0.75) for the inci- 
dent energy of 40 (15) AMeV. The process of evaporation 
(GEMINI) which affects the hot systems (SMF) is more 
relevant at higher incident energy. Thus, larger varia- 
tions of the isotopic composition (N/Z) of the largest 
residue (Ab; g ) as a function of the isotopic composition 
of the system manifest at 15 AMeV. Experimentally, 
FAZIA would measure such trend correctly. Concern- 
ing INDRA, the redundancy between the parametrisation 
used to deduce the mass number from the atomic num- 
ber (residue corridor), and the description of the decay 
by GEMINI globally results in an apparently good per- 
formance for stable systems; for such systems the evap- 
oration process ends in populating the residue corridor. 
However, such effect depends on the model. In particu- 
lar, for the most neutron-rich system, the mass number 
parametrisation reveals to be no more sufficient and a 
direct measurement of the masses is necessary [3^, [34{ . 
Such parametrisation will be even less adequate for exotic 
systems (which will be studied in future installations ex- 
ploiting exotic beams), when the residue corridor is not 
reached. To give an extreme example, we applied the 
experimental filters also to the hot systems (SMF), ne- 
glecting the decay. FAZIA measures correctly even these 
very neutron rich systems, while INDRA associates the 
residue-corridor masses to the atomic numbers and gives 
almost the same result which was previously simulated 
for the cold systems (SMF+GEMINI). 

In conclusion, in most studies of isospin effects on the 
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equation of state, or of thermodynamical properties in 
the decay, exotic systems are highly desirable; otherwise, 
the smearing effect of the evaporation process will be 
difficult to surmount. 



IV. THE EQUATION OF STATE: THREE 
EXPERIMENTAL APPROACHES 

The obscrvablcs studied in fig. [5] led to the conclu- 
sion that secondary decay may impose major difficulties 
in the experimental study of the equation of state. Ex- 
otic nuclei may be advocated as more suited systems, 
where the secondary decay may be better retraced as 
far as the residue corridor is not reached in the cooling 
process. Alternatively, we can also search for other ob- 
servables which are less sensitive to the secondary decay. 
Three possibilities are explored in this section: the im- 
balance ratio, first introduced by the MSU group [l7l.[l9j: 
the high-energy part of neutron and proton spectra, pro- 
posed in the framework of BUU [35(, and also success- 
fully exploited via IQMD [36| ; and finally, the neck emis- 
sion in ternary events, studied at length by the Catania 
group @,[33- 

A. Imbalance ratios 



However, in each rapidity bin, this observable is expected 
to be correlated to its value prior to decay. From these 
data points fig UJ is constructed. We may remark that, 
since the isotopic composition of the asymmetric system 
is intermediate between the symmetric systems, accord- 
ing to the definition of EqfJJ the imbalance ratio should 
approach zero in the vicinity of midrapidity and evolve 
towards negative values in the projectile side for increas- 
ing rapidity and eventually reach the value of —1 for 
large transparency; in the target side, the imbalance ra- 
tio should evolve for increasing negative rapidity towards 
the same values explored in the projectile side with op- 
posite sign. This overall behaviour, reflected by the hot 
fragments (SMF), is well reproduced also after the evapo- 
ration process (SMF+GEMINI); this shows that, indeed, 
the effects of secondary decay can be removed when suit- 
able obscrvablcs which combine differences and ratios of 
isospin dependent properties are adopted, such as the im- 
balance ratio. The difference between the predictions of 
the asy-stiff and asy-soft parametrisations is preserved 
after the secondary decay stage and is measurable with 
FAZIA. Comparing fig. [5] and fig. [7J we can see that the 
discrimination between the two parametrisations is es- 
sentially due to the presence of the neutron rich system 
68 Ni+ 68 Ni, which is the only case where the residue cor- 
ridor is not reached during the deexcitation. Following 
the discussion of fig. [5l we expect that the imbalance ra- 



Thc isospin-transport ratio, introduced by F.Rami [38[ 
to study isospin equilibration, has been exploited by the 
MSU group fiol [l7l [l9j because of its sensitivity to the 
density dependence of the symmetry energy. In the case 
of the Ni isotopes of the present study, the imbalance 
ratio is defined as 
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where x is an isospin sensitive observable, possibly lin- 
early correlated with asymmetry, measured in the three 
different reactions. The motivation of introducing such 
observable is that, by focusing on the differences in 
isospin observables between mixed and symmetric sys- 
tems, R is expected to largely remove the sensitivity to 
preequilibrium emission and enhance the sensitivity to 
isospin diffusion between projectile and target. We have 
undertaken this analysis for the reactions at 40 AMeV, 
fow which the statistics of binary events is good enough. 

Considering all impact parameters, we assigned the 
average isospin content of the largest fragment, not pro- 
duced in fusion events, to the x variable. This observable, 
explored in fig|3l is shown again in figlH together with 
the corresponding response of the experimental device 
FAZIA. From a closer inspection of figs l3!6[ one may no- 
tice that the effect of the secondary decay on this observ- 
able is not the same in the different rapidity bins. This 
can be ascribed to the fact that the excitation energy 
of primary fragments is not constant, but it is larger at 
mid-rapidities, which correspond to more central events. 



(N/Z) (A big ) 



1.16 



1.14 



1.12 



1.1 



1.08 



1.06 




-, • asy-stiff 
-, o asy-soft 



58 Ni+ 58 Ni 
40 AMeV 



-0.5 0.5 1 

K [relative units] 

z o.m. Ied 



FIG. 6: (Color online) Same observable studied in figO we 
focus on the cold systems (SMF+GEMINI) of the system at 
40 AMeV and we add the filtered data, simulated for the 
experimental device FAZIA. 
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FIG. 7: (Color online) Imbalance ratio (see text) for the 
asy-stiff and asy-soft parametrisations, corresponding to the 

same observable studied in fig[3] for the systems 58 Ni+ 58 Ni 
58 Ni+ 68 Ni and 68 Ni+ 68 Ni at 4Q A Me y ( for all i mpact param- 
eters, and as a function of the longitudinal component of the 
rapidity in the centre of mass relative to the projectile (tar- 
get) rapidity for positive (negative) rapidities. The statistical 
uncertainties of the model are indicated by the error bars. 



tio will provide a better discrimination between different 
effective interactions with the use of more exotic beams. 

A similar analysis could be performed taking the trans- 
verse energy E tT of light particles as a centrality selector 
(see fig. 2]) instead of the rapidity of the largest frag- 
ments. However, the latter appears more directly con- 
nected to the dissipation mechanism between the two 
reaction partners, while E tr may be influenced by ad- 
ditional effects, such as pre-equilibrium dynamics and 
cluster-emission mechanisms. For this reason, the analy- 
sis proposed in fig. [7] should be more helpful when com- 
paring experimental data with the predictions of trans- 
port models which differ by the treatment of correla- 
tions ,27}. Also different isospin-sensitive observables 
may be used, such as the isoscaling parameter and the 
isobaric yield ratio ln(y( 7 Li)/y( 7 Bc)) employed by the 
MSU group [nj, or the average N/Z of the light-charge- 
particle emission (l8j . This would be possible within our 
analysis protocol, but it demands a larger statistics for 
the secondary-decay treatment. For all those choices, as 
far as the observables are linearly correlated to the N/Z 
of projectile-like and target-like fragments [T3, EH, the 
corresponding results for the imbalance ratio should be 
rather close to the ones depicted in fig. [7] for the observ- 
able we adopted. 



B. Light particles: neutron and proton spectra 

The evolution of isotopic ratios of light particles as a 
function of their kinetic energy and the isotopic composi- 



sensitive to the stiffness of the equation of state [35|, |36| . 
When light particles (like t and 3 He) are concerned, the 
SMF model can no more access such observable which, 
in general, is strongly dependent on the model. For these 
reasons, in this section we restrict to the study of pree- 
quilibrium neutrons and protons. 

The neutron and proton yields as a function of the 
perpendicular component of the kinetic energy E±, are 
illustrated in fig. [g] for the system 58 Ni+ 68 Ni at 40 
AMeV. The spectra are simulated before (SMF) and af- 
ter (SMF+GEMIN1) the secondary decay. As expected, 
the low-energy side of the spectra is mostly affected by 
the secondary decay. To minimise the effect of secondary 
decay, we select the high-energy side of the spectra where 
the two simulations, before and after evaporation, coin- 
cide: we impose E± > 25MeV. 

The shape of the proton spectra can be measured di- 
rectly by an experimental device like FAZIA, as simu- 
lated in fig. [51 or like INDRA. The shape of the neu- 
tron spectra can be deduced from the proton spectra 
only if complete events can be recorded; since this is 
practically impossible to achieve at intermediate ener- 
gies due to the high multiplicity of the events, the neu- 
tron spectrum should be measured directly with the use 
of a neutron detector; in the present simulation, shown 
in fig. [51 we suppose that the detector DEMON is em- 
ployed. In particular, we suppose to dispose of DEMON 
telescopes covering all azimuthal angles, and only an in- 
complete interval of 30 deg in zenith angle, so that the 
shape of neutron spectra can be measured correctly. For 
such a configuration, we do not take into account the 
effect of the charged-particle detector on the neutron de- 
tection, even if, in a realistic configuration, such effect 
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FIG. 8: (Color online) Proton (left panel) and neutron (right 
panel) yields as a function of E± for the system 58 Ni-|- 68 Ni 
at 40 AMeV. Units are in percent of the total neutron yields. 
The spectra are simulated for the hot system (SMF - asy-soft) 
and for the cold system (SMF+GEMINI). The same spectra 
are simulated as measured with an experimental device com- 
posed of FAZIA coupled with DEMON modules. 
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could be large and the neutron detectors may have to be 
even substituted to some charged-particle detectors. It 
is possible, but it should be proved experimentally, that 
analogous information may be deduced from the analysis 
of £- 3 He spectra, which would make this analysis accessi- 
ble to neutron-blind devices like INDRA or FAZIA alone. 
The integral of the proton spectra can be measured di- 
rectly, by selecting only complete events. The integral of 
the neutron spectra can not be measured directly in any 
case: in general, the efficiency of neutron detectors would 
not be sufficient to impose a condition of event complete- 
ness even when using a 47r-array. The integral of neutron 
spectra is deduced by difference, from the measurement 
of the nuclear charge and mass of fragments and particles 
detected in complete events. Hence, for a correct evalu- 
ation of the efficiency of the neutron detector, also the 
mass detection, as provided by FAZIA, is necessary. The 
slight discrepancy between the theoretical and measured 
spectra in the simulation of fig. [HJ both for neutron and 
proton spectra, is due to the choice of defining as com- 
plete events those events where less than 15% of the total 
nuclear charge is missing. 

From the neutron and proton spectra we deduce the 
evolution of the ratio between the number of neutrons 
over the number of protons N„/N p , averaged over com- 
plete events and normalised to the isotopic composition 
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FIG. 9: (Color online) Ratio between the number of neutrons 
over the number of protons averaged over complete events 
and normalised to the isotopic composition of the system, as 
a function of the isotopic composition of the system. This ob- 
servable is simulated for incident energies of 15 and 40 AMeV, 
for hot systems, before the secondary decay occurs (SMF) and 
for cold systems, after the secondary decay (SMF+GEMINI). 
The two prescriptions for the equation of state, asy-stifr and 
asy-soft are employed. Specific event selections are imposed 
depending on the incident energy. 15 AMeV: only neutrons 
and protons with E± > 15 MeV are chosen. 40 AMeV: only 
neutrons and protons with E± > 25 MeV are chosen. In both 
cases, only events not leading to fusion are selected. 



of the system (N/Z) systenl , as a function of the isotopic 
composition of the system. In fig. |H] this observable is 
studied for hot systems, before and after the secondary 
decay, for the two incident energies 15 and 40 MeV per 
nucleon. We mention that the condition E± > 15 MeV 
is used for the reactions at 15 AMcV, due to the lower 
energy available in this case for the pre-equilibrium emis- 
sion. It has been already observed that the best sen- 
sitivity to the equation of state is obtained in central 
collisions J3j|. However, especially at 15 MeV per nu- 
cleon, it is difficult to extract central events from the 
E tI distribution, see fig. 2J Hence, for both incident 
energies, we consider all impact parameters, but we se- 
lect fusion events out. The different parametrisations 
for the equation of state, asy-stiff and asy-soft, do not 
determine any appreciable difference at 15 MeV per nu- 
cleon. This could be due to the fact that, at this low 
energy, compression-expansion effects are quite reduced 
and, while pre-equiquilibrium emission takes place, the 
nuclear density keeps rather close to the saturation value, 
where the two parametrisations give, by construction, the 
same value of the symmetry energy. On the other hand, 
at 40 MeV per nucleon, the two parametrisations of the 
symmetry energy lead to clearly distinct behaviours even 
for non-exotic systems: the normalised average neutron- 
to-proton ratio is up to 13% larger for the asy-soft case. 
Now the isotopic content of the pre-equilibrium emission 
is sensitive to the low-density behaviour of the symmetry 
energy, since particles mostly escape while the composite 
nuclear system is expanding. A neutron-richer emission 
is seen in the asy-soft case, corresponding to the higher 
value of the symmetry energy below normal density (sec 
fig-Hi)- One can also notice that the dependence of this 
observable on the system initial asymmetry is flatter for 
the reactions at 40 AMeV. This reflects the higher en- 
ergy available, with respect to the 15 AMeV case, at the 
pre-equilibrium stage. In the following, we concentrate 
on the systems at 40 MeV per nucleon. We observe that 
the signals are still preserved after the secondary decay, 
thanks to the energy selection we adopted. In fig. [TOl 
we simulated the response of the two detection devices, 
FAZIA and INDRA. The error bars account for the ac- 
curacy in the correction for the efficiency of the neutron 
detection, supposing that it is known with an uncertainty 
of around 5%; additional systematical uncertainties (not 
shown) would be produced by the procedure of account- 
ing for the number of fragments and particles which are 
not detected. 

The simulation demonstrates that FAZIA would per- 
form efficiently and succeed in measuring the signal. 
However, after the deexcitation stage, the mass of the 
evaporation residues, which should be measured in or- 
der to reconstruct the number of emitted neutrons, ap- 
proaches the residue corridor. In order to study a more 
constraining situation, for testing purposes, we replace 
the cold system by the hot system, before the decay, 
in the experimental filter. For the most neutron-rich 
system, we add simulated points representing INDRA, 
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FIG. 10: (Color online) Same observable as in fig. [5] for 
an incident energy of 40 AMeV. Both hot systems (SMF, 
left panel) and cold systems (SMF+GEMINI, right panel) 
are simulated as measured with an experimental device com- 
posed of FAZIA coupled with DEMON modules before and 
after secondary decay. The simulation of INDRA coupled 
with DEMON modules is added for the most neutron-rich 
system, before secondary decay. 



which docs not measure the masses, but deduces them 
through the evaporation-corridor prescription introduced 
in the detector simulation: these points, as they are re- 
ferred to a very neutron-rich system, far from the residue 
corridor, result in fact to be incompatible. We conclude 
that such signal survives after the secondary decay and is 
measurable only with a device like FAZIA in conditions 
where the reaction products are far from the residue cor- 
ridor. 



C. Neck emission in ternary events 

It has been discussed in the literature [|| that the 
dissipative dynamics leading to neck formation in mid- 
peripheral collisions at Fermi energies is strongly sensi- 
tive to the isospin transport properties and to the isovec- 
tor properties of the equation of state. 

We focus on the system 68 Ni+ 68 Ni at 40 AMeV; such 
incident energy is large enough to exceed the threshold 
for the formation of neck fragments. For this specific sim- 
ulation the stopping time was fixed to 160 fm/c, which is 
long enough to follow the dynamics of the neck formation. 
We select ternary events where, in addition to light par- 
ticles, we observe only the following three fragments: one 
quasi-projectile residue, one quasi-target residue, and one 
intermediate-mass fragment at midrapidity, which corre- 
sponds to the neck fragment. In addition, we impose 
two more constraints for the selection of ternary events. 
First, we reduce to semipcriphcral collisions where & re( j 
varies in the interval 0.45 < 5 rc d < 0.75: in fig. [5] we 



already signalled the presence of neck events in this in- 
terval. The second constraint is on the mass number of 
the third fragment: by adjusting a selection threshold 
to Z > 5, we ensure that the ternary events with three 
IMF's that we select can not be confused with secondary- 
decay events which produce fragments of small size. 

The observable we investigate is the isotopic composi- 
tion of the neck fragment in ternary events, selected as 
defined above. This system, studied before the secondary 
decay, is less neutron rich when the symmetry energy is 
parametrised with the asy-soft term, and more neutron 
rich when an asy-stiff interaction is used. This is due to 
isospin migration effects from projectile and target to- 
wards the low-density region of the neck. In particular, 
the evolution of the isospin asymmetry of the neck re- 
gion is conditioned by the contribution of the migration 
process to the currents of neutrons and protons, driven 
by density gradients. Such contribution is represented in 
cq. |31 by the term proportional to the derivative of the 
symmetry energy and the density gradient. The process 
of migration determines an enrichment in the neutron 
content of the diluted neck region, which is larger for 
larger variations of the symmetry energy as a function of 
the density, and is therefore more effective for a stiffcr 
form of the EOS around normal density , which is in this 
case given by the asy-stiff parametrisation (see fig. [IJ. 
When the secondary-decay process is added up, in the 
asy-soft case, the more excited neck fragments end up in 
the evaporation corridor and their N/Z ratio reduces. In 
the asy-stiff case, the hot system is more neutron-rich and 
the decay path does not reach the evaporation corridor. 

Such difference in shown in fig. 111! where the iso- 
topic composition of neck fragments in ternary events 
is plotted as a function of 0am, where </>am is the small- 
est angle which measures the angular misalignment be- 
tween the velocity vector of the neck fragment in the 
centre of mass frame and the direction along which the 
two heaviest fragments (representing approximately the 
quasi-projectile and the quasi-target) are aligned. This 
convention, followed in this paper, does not correspond 
exactly to the definition of p i a nc [2Q| , more often followed 
in experimental studies of the phenomenon of nonequi- 
librium fission. We signal that such studies are closely 
compatible with the present work and explicitly oriented 
to the investigation of the EOS. 

The increasing trend observed in the asy-stiff case may 
be due to the impact parameter mixing in our analysis. 
Indeed, fragments with larger transverse velocities are 
more easily emitted in the most central events of the se- 
lected interval of impact parameters, for which, due to 
the longer reaction time, isospin migration is more effec- 
tive. The trend is less clear in the asy-soft case, where 
isospin migration effects are less pronounced. For the 
same reason, fragments having 4>am close to 90 deg are 
more excited. In the asy-stiff case, the isotopic compo- 
sition of neck fragments reduces almost systematically 
when passing from the system before decay (SMF) and 
after decay (SMF+GEMINI). The effects related to the 
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different excitation energy are more evident in the asy- 
soft case: the less excited fragments experience a short 
decay path and are close to the isotopic composition of 
the hot fragments, while the fragments which could spend 
more excitation energy in the decay reached the evapo- 
ration corridor. The further addition of the experimental 
filter demonstrates that FAZIA would measure this ob- 
servable precisely. On the other hand, INDRA would not 
succeed in measuring this observable, and the discrep- 
ancy is evident when the residue corridor is not reached, 
that corresponds to the asy-stiff case. In conclusion, the 
asymmetry of the neck fragments appears as a promis- 
ing observable to be exploited in the investigation of the 
behaviour of the symmetry energy below saturation den- 
sity. Due to the isospin-migration mechanism, these frag- 
ments are particularly neutron-rich, especially in the asy- 
stiff case, and the effect survives to the secondary-decay 
stage. 

Isospin migration is more effective in neutron-rich sys- 
tems, however interesting results have been reported also 
in the case of stable, less asymmetric systems, in theo- 



b/j 

a 

,S-H 



o 
a 



1.5- 
1.4- 
1.3- 
1.2- 
1.1 
1- 



SMF 



asy-stiff 



asy-soft 



SMF + GEMINI 
INDRA 



SMF + GEMINI 



asy-stiff 



asy-soft 




SMF + GEMINI 

FAZIA 




20 ' 40 60 ' 80 



20 ' 40 60 80 



0AM 



FIG. 11: (Color online) Isotopic composition of the neck frag- 
ment in ternary events, selected as defined in the text, as a 
function of 0am, for the system 68 Ni+ 68 Ni at 40 AMeV. The 
spectra are simulated for the hot system (SMF, upper left 
panel) and for the cold system (SMF+GEMINI, upper right 
panel). The spectra representing the cold system are simu- 
lated as for the two experimental approaches INDRA (lower 
left panel) and FAZIA (lower right panel). 



retical [39(, as well as in experimental analyses [20(. It 
would be extremely interesting to pursue the investiga- 
tion of the neck dynamics also at lower beam energy. 
Though the analysis is complicated by the reduced statis- 
tics of ternary events, fragments are expected to be less 
excited, opening the possibility to access more directly 
the genuine dynamical effects. Work is in progress in 
this direction. 



V. CONCLUSIONS 

In summary, we built a protocol for simulating re- 
action experiments which profit from next-generation 
radioactive-beam facilities and from innovative experi- 
mental solutions; these latter are specifically dedicated 
to the measurement of isotopic observables in exclusive 
experiments at low and intermediate energies. In this 
work, the systems 68 Ni+ 68 Ni, 58 Ni+ 68 Ni and 58 Ni+ 58 Ni 
at 15 and 40 MeV per nucleoli were simulated for a 
continuous distribution of impact parameters b within 
the Stochastic Mean Field Model, with the addition of 
a secondary-decay process (described within the model 
GEMINI). The experimental conditions were simulated 
by supposing that the experimental observables are de- 
livered either by the former-generation detector INDRA, 
or by a next-generation detector, which corresponds to 
the FAZIA project; such project, under development, is 
intended to build a 47r-array of silicon detectors which, 
beside measuring the nuclear charge of the fragments, are 
also sensitive to the mass and have low thresholds. Such 
simulation protocol was used to discuss new experimen- 
tal solutions for probing the nuclear equation of state in 
heavy-ion collisions at Fermi energies. 

Firstly, we studied the effect of the secondary decay. 
Such process, which creates a significant bias in the de- 
termination of the EOS properties, should be taken into 
account in the choice of the experimental strategy. Ac- 
cording to our simulations, secondary-decay effects can 
be partially cured by studying exotic systems. However, 
even with exotic beams, secondary decay has to be con- 
trolled before a reliable extraction of the isovector equa- 
tion of state can be obtained. 

The strategy of constructing suitable combinations of 
isotopic observables (such as imbalance ratios) , that was 
proposed for similarly prepared systems, appears to be 
an efficient way to reduce the effects of the secondary 
decay. Other possible solutions consist in selecting ob- 
servables connected to fragments which are formed with 
insufficient excitation for feeding any relevant decay pro- 
cess, or uniquely connected to the initial phases of the 
reaction process, like in the preequilibrium stage. 

Following this second approach, we identified two more 
probes which are negligibly affected by the secondary 
decay: the relative yields of preequilibrium high-energy 
light particles (we focused on neutrons and protons) , and 
the isotopic composition of neck fragments in ternary 
events. Both observables would exhibit significantly dif- 
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ferent signatures for the asy-stiff or the asy-soft be- 
haviour of the symmetry energy, even when the sec- 
ondary decay is accounted for. This is particularly evi- 
dent in the case of the neutron-rich 68 Ni + 68 Ni system, 
for which isospin effects are enhanced. 

In this framework, observables related to the neck dy- 
namics appear rather promising. Interesting effects have 
been revealed already with stable beams (2Q, , which 
may constitute one passage towards exotic-beam experi- 
ments. Moreover, it would be very interesting to extend 
this analysis to lower beam energies. We found that these 
signatures could be identified with an exclusive detector 
device which is compatible with the FAZIA project, and 
which is based on the simultaneous detection of charges 
and masses over a large solid angle. The principle of such 
an experimental approach has several advantages. First 
of all, the possibility of performing a complete event re- 
construction through the data analysis allows to refine 



the criteria of centrality selection, which is a crucial point 
in the comparison between theoretical simulations and 
experimental data. Then, the measurement of the mass 
of all reaction products gives access to more sophisti- 
cated isospin-sensitive observables. Finally, we shall also 
stress that such experimental approach would be suited 
for analysing within the same experiment several probes 
for the symmetry energy, which are also sensitive to dif- 
ferent stages of the reaction. 
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